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We develop a non-perturbative numerical method to study a single electron tunneling through 
an Aharonov-Bohm ring in the presence of bound, interacting electrons. Inelastic processes and 
spin- flip scattering are properly taken into account. We show that electron-electron interactions 
described by the Hubbard Hamiltonian lead to strong dephasing and we obtain high transmission 
probability at $ = n even at small interaction strength. Depending on the many-electron state on 
the ring, dephasing can occur in elastic or inelastic channels with or without changing the spin of 
the scattering electron. 



Quantum interference can be studied in mesoscopic 
systems where the wave nature of electrons plays an im- 
portant role. Particularly noteworthy are studies of the 
Aharonov-Bohm (AB) oscillations in mesoscopic rings 
[El ||. Inelastic scattering of electrons is believed to be 
predominantly responsible for the loss of the phase co- 
herence in such experiments. When an electron interacts 
with optical phonons, the dephasing only occurs through 
inelastic processes |J. At low temperatures the phonon 
degrees of freedom freeze out, therefore other mecha- 
nisms for dephasing, like magnetic impurity scattering 
or zero-point fluctuations of the electromagnetic envi- 
ronment H have been proposed. It is nevertheless 
believed, that at low temperatures the electron-electron 
(e-e) interaction is a dominant mechanism for dephasing 

The AB geometries have been theoretically studied 
predominantly by self-consistent mean-field approxima- 
tions that break down for degenerate levels, which phys- 
ically happens at very low temperatures |t|, ||, ||. They 
do not describe transitions in which the symmetry of the 
wavefunction of bound electrons changes and they are 
thus inadequate to study decoherence. Renormalisation 
group techniques have been applied to AB systems [fTo) , 
where calculations were limited to spinless interacting 
quantum dot with two levels coupled to reservoirs. Par- 
ticular attention was devoted to the appearance of Kondo 
physics induced by changing magnetic flux, however due 
to limitations to spinless fermions, no spin-flip induced 
decoherence has been investigated by this method. 

To shed some new light on the problem of decoherence, 
there is obviously a demand for a capable method, that 
would treat the problem of the scattering of an electron 
through a limited region where e-e interactions would be 
exactly taken into account. In this Letter we propose a 
method that treats e-e interactions by direct diagonali- 
sation using iterative (Lanczos) technique. The method 
takes into account spin-flip processes, so it can also be 
used in calculations of spin-polarized transport JTT| . 

We apply the method to study a single-electron trans- 
mission through a ring with e-e interactions described by 
a Hubbard Hamiltonian. We show that dephasing can 



occur either by a) inelastic processes where the tunnel- 
ing electron excites bound electrons on the ring or by b) 
elastic processes, where the tunneling electron changes 
the symmetry of the degenerate many-electron wavefunc- 
tion. No exchange of energy is required in the latter case 
p] , |l2| , |l3| , p| : dephasing occurs because the tunneling 
electron leaves a trace on its "environment" , which con- 
sists of bound electrons. In either elastic or inelastic case 
the dephasing can occur with or without the spin-flip of 
the scattering electron. 

The proposed method is based on the multichannel 
scattering technique, that was developed for studying 
the tunneling of a single electron in the presence of 
scattering by phonons [|15| Since its introduction, 

it has been successfully applied to a variety of prob- 
lems where a single electron is coupled to phonon modes 
H [I?], |o[ |l| and even incorporated into Lan- 

dauer theory where the influence of electron-phonon scat- 
tering on the nonequilibrium electron distribution has 
been investigated We now generalise this method 

to study many-electron problems. While the method can 
be applied to more general situations, we choose for sim- 
plicity a particular physical system which will also serve 
as a toy-model for the study of the e-e interaction in- 
duced dephasing. We thus consider an AB ring coupled 
to two ideal one-dimensional leads, see Fig. (Q). The ring 
is described by a Hubbard-type Hamiltonian 
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The operator creates an electron with spin a at site 
j and we make a formal identification c-j^ = Ci t<7 K The 
phases <j)j describe phase changes due to magnetic flux 
penetrating the ring. The leads are described by the 
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Extended interaction region 




tron on the AB ring is given by 



Left lead Aharonov-Bohm ring 



Right lead 



FIG. 1: Aharonov-Bohm ring. Magnetic flux 3> penetrates 
the center of the ring. 



tight-binding Hamiltonian 

-fflead = — £lcad ^ (H+1,<T^ 0>i,<7 + h.C. 

-£icad E h+\,<x h,o + h.c. 

The operator creates an electron with spin a at site 
i on the left lead, while the operator bi^a does the same 
on the right lead. The ring is coupled to the electrodes 
with coupling constants t$: 



H c = -t ^2 {a\,<^ci,a + h.c.) - t ^2 [h 
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Our approximation consists of taking into account only 
those many-electron states in which at most one (scat- 
tering) electron is located outside the ring. Before the 
impact of the electron, there are n = nj + n± other elec- 
trons bound on the AB ring. We truncate all many-body 
states, where additional electrons hop from the interact- 
ing region to the lead. When physical parameters of the 
system, e.g. (e, t, U), are chosen in such a way that these 
n electrons are bound in the interacting region, the ap- 
proximation is equivalent to neglecting the exponentially 
decaying tails of the n-electron wavefunction in the leads. 
Before the scattering, the bound electrons are therefore 
in one of the n-particle eigenstates of the Hamiltonian 
i/ring, Eq. We denote these states by \a\) and their 
energies by ej. When an incoming electron with spin 
up enters the ring, the system is in one of the n + 1- 
particle states which we denote by |/3j). These states are 
not necessarily eigenstates of -faring- After the scattering 
there is a single electron in one of the leads, while the 
ring is again in one of the n-particle eigenstates of if ring- 
This state can either be one of the \a\) states or (if the 
spin of the scattering electron has been flipped) in one of 
the n-electron eigenstates with n-\ + 1 spin-up electrons 
and n± — 1 spin-down electrons. Spin-flipped states are 
labeled by \a\) and their energies by e\. 

By taking into account only the allowed states, the 
wave-function that describes the scattering of one elec- 
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We consider a steady-state scattering described by the 
the Schrodinger equation H\%>) = E\%>) with H = H lins + 
H\ cat i + H c . We operate on this equation from the left 
with (Pi | and using Eq. (|J) we obtain 

~*o E 6 L> d L> - *o E A> + E h ^k = Ee h 

j,cr j.cr k 

(4) 

where b's denote scalar products 6^- = (/3{|ci |0 -T|a?) 
and bf ja = (/3;|c 4i(J t |ap, while h Uk = (Pi\H ring \f3 k ) are 
the matrix elements of Hamiltonian H r { ng in the n + 1 
electron subspace. By operating with (aJ|ai )(T from the 
left we get 
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In an open outgoing channel (J, a) a plane wave can 
propagate, so that d\j = exp(ikj a )d\ • . By en- 
(2) ergy conservation the wave number kj a is obtained from 



e - 2ticad cos(K) = ej - 2ti ead cos(fej !(T ). The energy e 
is the initial energy of the n— electron bound state on 
the ring, K is the wave number of the incoming electron, 
and ej is the final energy of the bound electrons. Eq. (|5|) 
can thus be written as 



(E - (° + t load exp(ikj >a )) ' 



(6) 



Similar equation can be obtained for exponentially de- 
caying outgoing channels that we also take into ac- 
count. These are defined through the relation eo ~~ 
2t load cos(K) = cj - 2t lead cosh(fc J>0 .) and d\ j a = 
exp(— Kj,o)d\j a . Using this and Eq. ( |oj), t he leads can 
be removed (pruned) from the problem |l6| . 

At zero temperature, the electron scatters on the 
ground state of the n— particle state in the ring, |aj) with 
the energy eo. In the incoming channel we have both the 
incoming and outgoing waves, * = exp(— iKm) + 

rexp(iKm). We obtain (i^ot = ex-p(iK)d\ Q t + 
cxp(iK) — exp(— iK). The equation for the incoming 
channel thus contains an additional inhomogeneous term 
exp(iK) — exp(— iK). 

Equations ([|) , (||) and equivalent equations for the out- 
going channels in the right lead form a system of equa- 
tions for unknowns d\j >a) an d This system 
is solved for different energies of the incoming electron 
using the conjugated gradients squared (CGS) method 
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with a Jacobi preconditioner. The partial transmit- 
tivity through channel (J, a) is given by Tj. a (E) = 
sin(fcj )CT )/ sin(K)\dfj a \ 2 . Since the method is based on 
exact solution of many-electron problem, we can compute 
transmission at arbitrarily large values of U. 

Results can be improved by extending the interac- 
tion region, which is solved numerically by the Lanczos 
method, by adding additional sites from the leads. This 
procedure increases the computational Hilbcrt space, and 
consequently it properly takes into account decaying tails 
of bound electron wavefunctions in the leads. These im- 
provements mainly lead to energy shifts of the resonance 
peaks while the general characteristics of the spectra re- 
main unchanged. In principle, the region could be ex- 
tended until the desired convergence is achieved. In our 
calculations the interacting region consisted of the ring 
and one additional site from each lead. In cases where 
the ground state of the interaction region was degener- 
ate, we averaged the transmittivity spectra over all the 
degenerate states. The variational space taken into ac- 
count in our calculation was equivalent to a Hubbard 
model on 8 sites with no translational symmetry. The 
system of equations was solved to an accuracy finer than 
the linewidth in the calculated spectra. 

We now investigate the effect of the interactions on 
an electron as it tunnels through the ring. The on-site 
energies are e = — 4.5ti oa d, and the overlap integrals are 
t = \/3iicad, and we set £i ca d = 1- We have limited the 
energy of the incident electron to a half of the bandwidth, 
e.g. E = [— 2, 0], in order to avoid ionisation of the bound 
electron on the ring which would lead to two electrons 
exiting the interaction region. In all cases, the incoming 
electron was chosen to have spin up. 

For a test case we first consider a single electron with 
spin down bound on the ring. We start with the nonin- 
teracting case. In the absence of the magnetic field the 
transmission reaches unity at the resonance, Fig. (|2]a,) . 
The electron is fully reflected at any incident energy when 
the magnetic flux is $ = n, Fig. ( ||b). This is the usual 
Aharonov-Bohm effect. 

We now turn on the interaction. At <I> = we 
still see a unitary peak at the energy of the single- 
electron resonance, followed by smaller satellite peaks 
caused by the interaction, Fig. (||c). At $ = w, 
when in the absence of U the electron is fully re- 
flected, we obtain very high transmission probability 
despite relatively small U — 0.2, Fig. (||d). In the 
largest peak the transmission approaches the value 1/2. 
Since the incoming electron and the bound electron 
are not entangled, their total spin is not well defined, 
therefore the total wavefunction is a superposition of 
a singlet and a triplet state with equal amplitudes: 
|U) = l/V2(\S=l,S z = 0) + \S = 0,S z = 0)). The 
triplet scattering has zero transmission probability at 
$ = 7r since two electrons with spin-up do not inter- 
act. The singlet scattering, however, reaches the unitary 
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FIG. 2: Transmission probability as a function of the incident 
electron energy for one electron with spin down bound on the 
ring. The incoming electron had the spin up. The coupling to 
the lead is to = 0.4. In all cases transmission is purely elastic. 



limit at finite U at the main resonance peak. Averaging 
over both contributions, we get T — 1/2. 

The spin-flip scattering part of the transmission proba- 
bility is shown in the inset in Fig. (^d) . The spin- flip and 
normal scattering contribute equally to the total trans- 
mission probability. Both are purely elastic with respect 
to energy changes. 

Transmission in the spin-flip channel occurs because a 
spin- flipped electron loses its phase memory Jll[]. Trans- 
mission without the spin-flip occurs because of the dou- 
ble degeneracy of the ground state of the bound elec- 
tron. The tunneling electron changes the symmetry of 
the bound electron wave- function and it thereby acquires 
a phase shift. Such symmetry-changing transitions are 
only possible when the electrons interact. 

We now focus on the case of several strongly interacting 
(U = 1) bound electrons when the flux is $ = ir. When 
the bound state consists of three electrons with spin up 
and one electron with spin down (nj = 3,n^ = 1, see 
Fig. (||)) no spin-flip scattering is possible because such 
processes are energetically impossible. The ground state 
is however fourfold degenerate and the tunneling electron 
can get through the ring at finite U by changing the sym- 
metry of the many-electron state on the ring. Since the 
ground state is degenerate, this process is purely elastic. 
In the case of n-f = 0, rij, = 4 the ground state is non- 
degenerate, however the spin-flip processes are energet- 
ically allowed. We therefore obtain transmission proba- 
bility only in spin-flipped channels. Since in this case the 
ground state is not degenerate, the transmission consists 
of purely inelastic processes. 

In the case when the ground state is degenerate and 
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case of n 
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2, where at U = 1 due to widely 
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FIG. 3: Transmission probability as a function of the incident 
electron energy for 71 f (nj,) electrons with spin up (down). 
Incoming electron had spin up. Interaction U = 1.0, coupling 
constant to = 0.3. 




spaced many-electron levels transmission remained zero 
in the whole interval of incoming electron energy. At 
large U — 15 the energy difference between the nondegen- 
erate ground state and the first excited state decreases 
in comparison with [7 = 1 case, as the states become 
compressed in the lower Hubbard band. We changed the 
on-site energy to e = —20 in order to keep the electrons 
bound on the ring. At $ = there are several energies 
at which the electron can resonantly tunnel through the 
ring, Fig. (^). At $ = n, the electron can only tunnel 
inelastically. The energy difference to the first excited 
state in the n electron Hubbard band is approximately 
1.4. We find indeed that only the electrons that are more 
than 1.4 above the bottom of the energy band can tunnel, 
Fig. (0b). Such inelastic processes occur both without 
(Fig. Qc) or with spin- flip (Fig. ||d). 

Using a simple model and a new numerical method we 
have extracted the essential physics of tunneling through 
the AB ring. In particular we have focused on the the 
role of e-e interactions on dephasing. While the proposed 
method clearly has some limitations (small interacting re- 
gions, inability to describe ionization processes, neglect 
of many-body effects in the leads), it nevertheless gives 
precise answers to the question: 'What are the main de- 
phasing mechanisms caused by the e-e interaction?' 

A particle can tunnel through AB ring at $ = n elasti- 
cally by a) changing the symmetry of the many-electron 
state which is possible in the case of degeneracy or b) by 
flipping the spin. Tunneling can also occur in the inelas- 
tic channel by exciting the many-electron state on the 
ring into an excited state with or without the spin-flip. 
Depending on the number of bound electrons, their total 
spin, degeneracy of the ground state and available energy 
of the incoming electron, the total transmission can be 
composed of partial transmissions caused by either one 
of the listed processes. 

Authors acknowledge the support of the Ministry of 
Education, Science and Sport of Slovenia. 
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FIG. 4: Transmission probability as a function of incident 
electron energy, for rif = 2, nj, = 2, U = 15 and e = —20. 

the spin-flips are allowed, we expect dephasing to occur 
with or without spin flip. Such is the case of nj = 1, n± = 
3. The transmittivity without spin-flip is purely elastic, 
while the spin- flip processes are predominantly elastic, 
with small contribution from inelastic channels. Finally, 
for n f = 2,n± = 2 electrons are fully reflected from the 
ring since there are no allowed spin-flip nor symmetry- 
changed channels. 

Finally we show the influence of large U = 15 on the 
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